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Recent results on QCD studies in e + e~ annihilations at LEP are presented. Data recorded by the LEP 
experiments at centre-of-mass energies between 91.2 to 206 GeV are included. The main topic is the measurement 
of a s from event shape variables and associated aspects like the energy evolution or non-perturbative power-law 
corrections. These 'standard' measurements are complemented by new determinations using the 4-jet rate with 
an excellent precision. A summary of the results on QCD colour factors from angular correlations in the 4-jet 
system completes this report. 



1. Introduction 

The LEP experiments ALEPH, DELPHI, L3 
and OPAL have collected 700 pb~l of annihila- 
tion data from 1989 to 2000. Final QCD analyses 
of these data using the latest detector simulation 
and advanced correction techniques are now being 
published. The experiments presented cumula- 
tive summary papers including various centre-of- 
mass energies and a wide collection of observable. 
The experimental systematic uncertainties are for 
many observables as small as 1%, and the statis- 
tical errors typically 0.1 %. Hence, these data of 
un-preceeded precision will serve as reference for 
future experiments. This report summarises re- 
cent measurements of a s from event shapes and 
4-jet observables, investigations of power law cor- 
rections and results on QCD colour factors. 

2. Measurements of a s from event shapes 

Theoretical calculations for event-shape dis- 
tributions are available at next-to-leading or- 
der (NLO) complemented by all-orders resumma- 
tion of large leading and sub-leading logarithms 
(NLLA) for certain observables. A unified pre- 
scription for the matching of fixed-order and re- 
summed calculations, the so-called modified logR 
matching scheme, has recently been suggested L 
and is applied by the LEP experiments for their 
analysis. The most commonly used variables are 
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thrust, heavy jet mass, wide and total jet broad- 
cnings, C-parameter and j/3, the 3-jet resolution 
parameter in Durham scheme. These observ- 
ables were also selected by the LEP QCD work- 
ing group for a LEP combination . A virtue of 
event-shape distributions compared to other ob- 
servables is that even with limited event statistics 
at LEP2 energies a measurement of a s is possible 
with a reasonable precision and enables the obser- 
vation of the energy evolution of a s . An exam- 
ple of a final analysis is shown in fig. ^ where the 
measurements of the wide jet broadening variable 
at all LEP energies are compared to the result 
of fits with the NLO+NLLA theoretical predic- 
tion. ALEPH has in addition applied a pertur- 
bative NLO correction for the b-quark mass. In 
general the description of the data is good, al- 
though restricted to the central part of the dis- 
tributions at LEP1 where the high precision of 
the data requires NNLO calculations. The mea- 
surements have been corrected for acceptance and 
detector resolution effects and the fits are carried 
out at hadron level. Therfore, the perturbative 
prediction is folded to the hadron level by means 
of transition matrix accounting for the hadro- 
nisation, obtained from standard QCD genera- 
tors like PYTHIA, HERWIG or ARIADNE. The 
systematic uncertainties for the measurement of 
a s is dominated by theoretical uncertainties in- 
duced by missing higher orders. The LEP QCD 
working group advises the recommendation of 
for the estimation of perturbative uncertainties. 
This method combines different estimates in the 
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Figure 1. Distributions at various centre-of-mass 
energies of the wide jet broadening compared to 
theoretical predictions used to determine a s . 



'uncertainty band method', which takes not only 
standard renormalisation scale variations but also 
a variation of xl , the resummed logarithmic vari- 
able re-scaling factor, into account 

2.1. Combined measurements 

The measurements using different variables and 
different energies are combined in single numbers 
per Ecu and finally using the QCD-predicted 
evolution in global result for a s (Mz)- In fig.Elthe 
combined result from L3 for a s (Q) is shown, in- 
cluding also measurements using radiative events 
resulting in reduced centre-of-mass energies be- 
low Mz- The experiments have applied different 
techniques for their combinations in terms of the 
correlation of systematic uncertainties, but in all 
cases the theoretical uncertainties appear to be 
largely correlated, both between different vari- 
ables and between different energies. As a con- 
sequence, the gain in precision of the combined 
measurements in limited and the combined un- 
certainty appears as average rather than an im- 
proved uncertainty. 
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Figure 2. Measurements of a s , combining five 
event-shape variables, as function of centre-of- 
mass energy. 



ALEPH has carried out a global analysis of 6 
variables at 8 energies between 91 and 206 GeV 
0], based on an integrated luminosity of about 
700 pb _1 . Their combined result for a s (M z ) is a 
weighted average with weights proportional to the 
inverse square of the total individual errors. To 
account for correlations of systematic uncertain- 
ties the whole combination is repeated separately 
for all the individual variations of the analysis. 
The final result is 

a s (M z ) = 0.1214 ± 0.00U exp ± 0.004%. 

L3 presented in a collective preprint 5 various 
studies of QCD, including measurements of a e 
from an earlier publication 6 . They combined 
the same five variables, but excluded —In 2/3. 
Measurements at reduced centre-of-mass energies 
using radiative events were included. Theoreti- 
cal uncertainties are estimated by a variation of 
the renormalisation scale and different matching 
schemes. L3 has chosen at LEP2 very large fit 
ranges, in contrast to all other experiments, in 
order to reduce the statistical uncertainty. This 
leads on the other side to larger theoretical er- 
rors. The combination method proceeds in two 
steps: first an unweighted average of the five vari- 
ables is build at each Eqm > second these combined 
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measurements are fit by the QCD evolution for 
a s (Mz). At that stage correlations are included 
by the assumption of minimal overlap of system- 
atic uncertainties. The result of this combination 
is: 

a s (M z ) = 0.1227 ± 0.0012 e:rp ± 0.0058^. 

DELPHI updated a publication *7, on event- 
shapes with a new re-analysis [S] according to 
the LEP QCD scheme. Traditionally DELPHI 
presents measurements based on three perturba- 
tive methods: the standard modified logR match- 
ing scheme, the pure NLLA calculation (valid 
only in a narrow 2-jet region) and the fixed or- 
der prediction alone with an optimised scale i° ?f . 
The perturbative uncertainties are estimated by a 
variation of xl only for the logR/NLLA schemes 
and by a renormalisation scale variation of a fac- 
tor of two around the optimum £° p *. The same 
set of variables as used by L3 is combined at en- 
ergies between 89 and 207 GeV, hence including 
also off-peak data at 89 and 93 GeV separately. 
The combination technique applied by DELPHI 
follows closely the LEP QCD Ansatz, but makes 
in addition the minimum overlap assumption for 
hadronisation and perturbative systematic uncer- 
tainties. The final result, split into the three per- 
turbative schemes, reads as: 

LogR a s = 0.1205 ± 0.0020 ea: p ± 0.0050 t/l , 
0{a 2 s ) a s = 0.1157 ± 0.0018 ea; p ± 0.0027 t/l , 
NLLA a s = 0.1093 ± 0.0023 exp ± 0.0051*/, . 

OPAL is currently preparing a new analysis 
on a s from event shapes and preliminary results 
were already included in the LEP average. In the 
last OPAL publication 9 data from LEP up to 
189 GeV were combined with lower-energy data 
from JADE at 35 and 44 GeV. The variables are 
the differential 2-jet rate (equivalent to ys) and 
the mean jet multiplicity using the Durham and 
Cambridge jet-finding algorithms. The combined 
result is: 

a s (M z ) = 0.1287±0.0012 &Ep +0.0034 t , l -0.0016 t/l , 

where a significantly asymmetric uncertainty is 
observed for the scale variation. 



3. power law corrections 

Non-perturbative effects in hadronic observ- 
ables in e + e" annihilation are scaling with powers 
of 1/Q and can be described by analytical models 
of power law corrections |10j. Power corrections 
in the spirit of these models are related to in- 
frared divergences of the perturbative expansion 
at low scales. Analytical calculations introduce 
one additional phcnomcnological parameter ao, 

1 f* 1 

ao(w) = — a s (k)dk , 
Mi Jo 

which measures effectively the strength of the 
coupling up to an infrared matching scale /Lti of 
the order of a few GeV. The parameter ao is ex- 
pected to be universal and must be determined 
by experiment, usually in conjunction with a s - 

An improved theoretical prediction is obtained 
by merging perturbative and non-perturbative 
terms. This yields for event-shape mean values, 

(V) = (Vpert) + (2/powor), 

for a generic variable y, where the additive power 
correction term is given by 

(y P owor) = c y ■ P(ao)/Q, 

with a variable-dependent constant c y . In event- 
shape distributions the power correction appears 
as a shift of the perturbative spectrum by the 
same additive term 

D v (y) = D pett (y- C -^l). 

For the jet broadenings this shift is not constant 
but depends on the value of the broadening. In 
fig. |21 different mean values are shown as func- 
tion of y/s and compared to pure perturbative 
and power-corrected predictions. 

A good description of the data is achieved both 
for mean values and the central part of distribu- 
tions. Different groups have analysed often simi- 
lar datasets including lower energy measurements 
and determined the two parameters a a {M%) and 
ao(2 GeV) from a simultaneous fit. The com- 
bined final results are given in table The 
value of «o(2 GeV) is close to 0.5 and the value 
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Figure 3. Measurements of event-shape mean val- 
ues as function of centre-of-mass energy compared 
to fits including power corrections. 



of a s (Mz) around 0.114, significantly lower than 
with standard Monte Carlo corrections for hadro- 
nisation. The two parameters in the simultane- 
ous fit are strongly correlated. It is instructive 
to consider the results in the plane a s (Mz) ver- 
sus ao(2 GeV), shown for the ALEPH measure- 
ments in fig. 0] It appears that the jet broad- 
ening variables prefer a significantly higher value 
of cto and a lower value of a s , incompatible with 
thrust and C-parameter, and with a s determined 
with Monte Carlo corrections. The large system- 
atic uncertainty for ao from the jet broadenings is 
traced back to uncertainties of the perturbatively 
calculated power correction term, which might 



Table 1 

Combined results on a 3 (Mz-) and ao(2 GeV) using 
different variables. Partly mean values or distribu- 
tions were analysed. 

a s (M z ) a {2 GeV) 
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0.1112 
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0.0053 


0.496 


± 


0.101 


distr. 
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Figure 4. Contours of confidence level for simulta- 
neous measurements of a s and ao (ellipses) com- 
pared to the combined measurement of a s using 
Monte Carlo corrections (shaded band). 
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be responsible for the large spread between dif- 
ferent variables, spoiling the predicted universal- 
ity of ao- Furthermore the findings of different 
groups are not always consistent with each other, 
as demonstrated in fig. |U where central values 
and contour ellipses of ALEPH and Movilla et al. 
are consistent while the value of ao from DELPHI 
is much lower and the error ellipse significantly 
smaller. 
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Figure 5. Comparison of confidence level contours 
from different groups for determinations of a s and 
ao using the total jet broadening. 



4. 4-jet observables 

Traditionally angular 4-jet observables have 
been used to determine QCD color factors, often 
in conjunction with a determination of a s from 
the 4-jet rate. The precision of the measurements 



is substantially improved with the advent of NLO 
0(a>l) calculations. Two analyses have been car- 
ried out recently by the ALEPH [TT] and OPAL 
|12) collaborations, which determined simultane- 
ously the colour factor ratios Tr/Cf, Ca/Cf and 
the coupling constant. Using the QCD value of 
normalisation Tr=1/2, ALEPH obtained: 

C A = 2.93±0.14 stat ±0.58 sys , 
C F = 1.35± 0.07 5tat ±0.26 SHS , 
a s (M z ) = 0.119 ± 0.006 stat ± 0.026 sas . 

These measurements of the colour factors are 
excellent agreement with the QCD expectations 
C A = 3 and C F = 4/3. The OPAL analysis yields 

C A = 3.02 ± 0.25 stai ± 0A9 sys , 
C F = 1.34±0.13 itat ±0.22 s „ s , 
a s (M z ) = 0.120 ± 0.011 stat ± 0.020 sys , 

again in good agreement with the ALEPH mea- 
surement. The confidence level contours of the 
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tios compared to various gauge group values. 
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ALEPH and OPAL analyses are shown in fig. 
in the plane Tr/Cf versus Ca/Cf- The measure- 
ments confirm the QCD expectation of SU(3). 

4.1. a s from the 4-jet rate 

Alternatively, setting the colour factors to 
QCD values, a s can be determined from certain 
4-jet observables like the 4-jet fraction at NLO. 
A better sensitivity is obtained compared to 3-jet 
observables, since the leading term is already in 
a 2 s . For the 4-jet rate defined in the Durham and 
Cambridge schemes also resummed calculations 
are available, although only for the R matching 
scheme. ALEPH [TT] and DELPHI [T3] have pre- 
sented measurements of a s using this technique, 
a new OPAL analysis including also data at LEP2 
and from JADE is under way. 



inside the fit range provided that the renormal- 
isation scale is set the experimentally optimised 
value £° pt = 0.73. The quality of the fit becomes 
rapidly worse for scales off the optimum. Sys- 
tematic uncertainties for the 4-jet analysis are not 
dominated by perturbative but hadronisation un- 
certainties. This is related to the fact that 4- 
and 5-jet production are less well described in 
standard Monte Carlo generators. ALEPH ap- 
plied a Bayesian method to determine the size of 
systematic uncertainties, which consists of a de- 
weighting of models or theories yielding a bad de- 
scription of the data (large x 2 ) by re-scaling the 
resulting uncertainties. This procedures leads to 
very small errors. 



ALEPH 




Figure 7. The 4-jet fraction at detector level 
as function of lny cut compared to a fit of the 
0(c^)+NLLA prediction. 



ALEPH measured the 4-jet rate in the Durham 
at LEP1 only and determined a s from a fit of the 
C(a^)+NLLA (R matching scheme) prediction, 
as shown in fig. [7| The data are well described 
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Figure 8. Fits to the 4-jet rate measured at LEP1 
for two different scale evaluation methods. 
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DELPHI investigated the two jet finding 
schemes Durham and Cambridge and concluded 
that the Cambridge scheme has smaller system- 
atic uncertainties. Their measurement and fit 
result are shown in fig. |SJ where also the ra- 
tio of data over theory is given for = 1 and 
cc° pt . DELPHI used in contrast to ALEPH only 
the fixed order calculation and found a need for 
rather small scales, x° pt = 0.015 for Durham and 
£° pt = 0.015 for Cambridge. The perturbative 
systematic uncertainty was estimated by a vari- 
ation of a factor of 2 around x° pt , this gives rise 
to very small uncertainties. DELPHI confirms 
that the main systematic uncertainty is associ- 
ated with the description of R4 in the QCD gener- 
ators. The 4-jet measurements are summarised in 
table |21 the ALEPH result is for a direct compar- 
ison with DELPHI also given in a non-Bayesian 
approach for the systematic error. The total un- 



Table 2 

Measurements of a a from the 4-jet fraction. 





DELPHI 


ALEPH 


ALEPH 




Cambridge 


Durham 


non-Bayesian 




0.1175 


0.1170 


0.1170 


exp. 


0.0009 


0.0008 


0.0008 


had. 


0.0027 


0.0004 


0.0021 


pert. 


0.0007 


0.0009 


0.0016 


tot. 


0.0030 


0.0013 


0.0027 



certainties of about 3 % for these measurements 
are very competitive compared to methods using 
3-jct observables, in particular the perturbative 
uncertainties are at the 1 % level, similar to full 
NNLO determinations [T5| . 

5. Conclusion 

A wealth of measurements of hadronic observ- 
ables has been provided in 11 years of data taking 
at the LEP collider. These measurements allowed 
the LEP collaborations to perform detailed tests 
of perturbative QCD and determinations of the 
fundamental parameters. A collection of the most 
important results is shown in fig. [5| All collab- 
orations have provided measurements of a s from 
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Figure 9. Summary of measurements of a s (JV/z) 
at LEP using different methods and observables. 



event-shape distributions and combined individ- 
ual results from various variables at energies from 
Mz to 209 GeV. The systematic uncertainties are 
dominated by missing higher orders and in order 
to match the experimental accuracy of 1 % NNLO 
calculations are required. Event-shapes have fur- 
ther been investigated in the context of power-law 
corrections. Simultaneous fits have been carried 
out to determine the non-perturbative parameter 
ao and a s . While the qualitative description of 
the data is good, the quantitative interpretation 
is not unambiguous since the value of a s with 
power corrections is significantly lower than with 
Monte Carlo corrections for hadronisation. Fur- 
thermore, the predicted universality of cvo is ver- 
ified at the level of 20% of its precision only. In 
particular the jet broadening variables need to be 
further investigated. 
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Complete NLO calculations for 4-jet produc- 
tion improved significantly the determinations of 
the QCD colour factors, which are found to be 
in good agreement with SU(3). New measure- 
ments of a s from the 4-jet rate using these calcu- 
lations have been presented. The precision of this 
method is very good and yields a 1 % uncertainty 
for missing higher orders. The total uncertainty 
is of 3 %, dominated by model uncertainties. The 
4-jet method is comparable in precision with de- 
terminations using the fully inclusive observables 
Rz and R T , for which complete NNLO calcula- 
tions are available. 
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